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NEW EXAMPLES OF DETERMINANT DIVISIBILITY 

SEQUENCES 

KRZYSZTOF GORNISIEWICZ 


Abstract. In this paper we consider divisibility sequences obtained from square 
matrices. We work with of matrix divisibility sequences associated to a semigroup 
and arising from endomorphisms of an affine space. We prove that determinant 
divisibility sequences originated from powers of square matrices are generalized 
Lucas sequences. 


1. Introduction 

By the divisibility sequence we mean in this paper a sequence {d n }n e N of inte¬ 
gers such that if n\m then d n \d m . One of the most famous divisibility sequence is 
the Fibonacci sequence: 0, 1, 1, 2, 3, 5, 8, 13, 21, 34,... which arise from linear recur¬ 
rence: F n = F n _i + F n _ 2 . This is an example of the Lucas sequences: L n = a a Zp ? 
where a, /3 are the roots of some quadratic polynomial over Z. See [2] for a complete 
classification of linear recurrence divisibility sequences and [5], [6] for introduction 
to other divisibility sequences. In this paper we discus properties of certain matrix 
divisibility sequences. We follow the approach initiated in [1], 


2. Matrix divisibility sequence 

Let S' be a commutative ring with 1. Let M r (S ) be a ring of r x r matrices with 
entries in S. By a divisor class of a matrix M G M r (S ) we mean a coset GL r (S)-M of 
M with respect to the natural left action of GL r (S) . We say that matrix M G M r (S) 
divides a matrix N G M r (S) if there exists a Q G M r (S ) such that N = QM. If M 
divides N, then any element of the divisor class of M also divides N. Let (r, •) denote 
a semigroup. A divisibility sequence of matrices over a commutative ring S, indexed 
by T, is a collection of matrices {M a } aer in M r (S), such that if a divides /3 in T, 
then M a divides Mp in M r (S). If {M a } ae p is a divisibility sequence of matrices, then 
by the multiplicativity of the determinant {det(M a )} aeT is a divisibility sequence of 
elements of the ring S. 
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We fix a faithful representation: 

[■] : r -)• End(A™) : a -A [a] 

of T into the group of endomorphisms of affine m-dimensional space A™ over S. 


Definition 2.1. Let x e A r s . The matrix divisibility sequence associated to (T, [•]) 
is the sequence of Jacobians {J ct (a;)} Q , e r which are r x r matrices with (i,j)— entry- 
given by partial differentials: 

[Ja( x )ki ■= d(([a](x))i)/dx h 

where ([a](a;))j is an ith entry of the value of the endomorphism [a] on x. The asso¬ 
ciated determinant divisibility sequence is defined by {det(J a (x)} ae p. 


3. Main Result 

Theorem 3.1. Let X e GL r ( Z) and Ai,..., A r be eigenvalues of X. Then for every 
n > 1: 

( \n _ \n \ 2 

t— r 
A ' •' ' 

is an integer and the sequence { D n } ne pj is a determinant divisibility sequence. 


Proof: Let A", Y. Z be square s x s matrices. Assume that entries of matrices Y and 
Z are functions of entries of the matrix X. Then the following matrix derivative 
formula holds ([4]): 


(3.2) 


d{YZ) 

dX 


(/ <8> Y) 


dZ 

lx 


+ {Z f <8) /) 


dY 

dX' 


where <8) means the Kronecker product, / is the identity matrix of rank s and A* 
means the transpose matrix of A. In addition we will use property of the Kronecker 
product: 


(3.3) 


(.4 ® C)(B ® D) = AB ® CD 


for any square matrices A, B , C, D of size sx s. From now on we fix T = N. Consider 
the group G of all invertible s x s matrices with the embedding: 


G —)• A s : 


A n ••• A ls 
A.s i ■ ■ ■ X ss 


l— t (An, • • • > Ads, • • •, A s i,..., A S s) 
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We define the endomorphism [n] for n G N. Let X := [X tJ ] G G and respectively 
X n \= [Xki\ £ G, where we treat Xu as functions of A^-, for 1 < < s. We 

dehne [n] : A s2 —>■ A s2 as 

[u] (An, ■ ■ ■ j X \ SJ ..., X s i, • • •, X ss ) (Aii,... , Ai s ,..., A s i,..., A ss ) . 

Using (13. 2 p we compute Jacobians of the n-th power of the matrix X 


_ d(X n ) 
~ dX 


d{ X n - x X) 
dX 


(/ <g> A" n_1 ) 


dX 

dX 


+ (X* <8) I) 


dX n ~ l 

dX 


By induction and the property (13.3ft of the Kronecker product we get: 


n —1 

J n = Y^{X f ) k ® X 11 - 1 ^. 
k =0 

Assume that X G GL s {7j) is a diagonalizablc matrix . Then X = PDP -1 , for 

some P G GL S (C) and a diagonal matrix D = diag{X±, X 2 , ..., A s }. The set of 

matrices with distinct eigenvalues is dense in the set of all square matrices M s ( C) 

2 

with respect to the topology of C s . Hence we can assume that X has different 
eigenvalues. Therefore 


n—1 


j n = ^((PDP-'yy <g> (pdp 


— l\n—l—k 


n—1 


= > (p- i y(D k y(py®PD 


~\n—l—k t~)—1 


k =0 


k =0 


p~ L = 


= {(p- 1 y®p) 


n—1 


P k ® D n ~ x 


— 1 — k\ 


,k =0 


(P^P 


-1\ 


and the determinant divisibility sequence is of the form: 


D, n = det = det 


n—1 




,k =0 


n—l—k • 


Since D k ® D 


we conclude that: 


is a diagonal matrix whose diagonal consists of terms X- X 


k\n—l—k 


* 3 


n—1 


D„=n 2 n^r 1 n e A br i_ ‘ =" 2 \daxf- 1 n 

k =0 l<i<j<s 


1=0 


X i X 1 


The last term on the right hand side is a product of values of symmetric polynomials 
computed at eigenvalues of the matrix X. The Galois group of the splitting held of 
the characteristic polynomial of X acts trivially on these algebraic integers, hence 
D n G Z. For any n,m such that n divides m we have [det X ] n \[det X] m and (Af — 
A”) | (A™ — A" 1 ) . Therefore, the sequence D n is a divisibility sequence. For equal 
eigenvalues we compute D n using the exact form of symmetric polynomials instead 
of their fractional expression. 

If X is not a diagonalizablc matrix, then instead of the matrix D we consider 
an upper-triangular matrix obtained from the Jordan form of X with eigenvalues 
Ai, A 2 , • • •, A r on the diagonal and eventually we get the same sequence D n . 
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4. Examples 


1) Let X G GL 2 (Z) and a = trX, b = tr 2 A" — 4detAL Then using Theorem 1 we 
obtain the sequence presented in [1], example 4.3: 



2) Let X G GL 3 (Z) and b = — trX, c = Adi + X 22 + A 33 , d = — det Ad The 
discriminant of the characteristic polynomial of X is A = (4 Aq — A 2 )/27, where 
A o = b 2 — 3c and Ai = 2 b 3 — 9 be + 27 d. We obtain the divisibility sequence defined 
by: 


D 


n 


n 2 d n 1 

A 



b + e l A + e 2 l A 0 A\ 

3 ) 


^b + e i+1 A + e 2 * + 2 ) A 0 A^ n ’ 


where A = yj (Ax + y/— 27A)/2, A = (Ai — 1 /—27A)/2 and e is a fixed primitive 
cube root of unity. 

3) It is easy to compute values of D n for any square matrix X. The matrix X = 


1 -2 -6 
0 1 3 

-1 0 1 


G Gl 3 (Z) gives the following divisibility sequence: 


n 

dn 

factorization of d n 

1 

1 

1 

2 

100 

2 2 5 2 

3 

6561 

3 8 

4 

193600 

2 6 5 2 11 2 

5 

808201 

29 2 31 2 

6 

189612900 

2 2 3 8 5 2 17 2 

7 

50131657801 

41 2 43 2 127 2 

8 

4096576000000 

2 12 5 6 11 2 23 2 

9 

159625511221401 

3 14 53 2 109 2 

10 

1865976489302500 

2 2 5 4 29 2 31 6 

11 

31583922467632921 

131 2 857 2 1583 2 

12 

21985833099924302400 

2 6 3 8 5 2 11 2 17 2 71 2 109 2 

13 

2370466451421685365841 

1637 2 4057 2 7331 2 

14 

118070682478980566428900 

2 2 5 2 41 2 43 6 83 2 127 2 

15 

2362255369723766871090801 

3 8 29 2 31 2 2969 2 7109 2 

16 

84956038709284864000000 

2 1s 5 6 11 2 23 2 47 2 383 2 
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-12 4-1 

0 1-2 2 

4) The matrix X — ^ ^ e Gl±\L) gives the following divisibility 

0 10 1 

sequence: 


n 

d n 

factorization of d n 

1 

i 

1 

2 

65536 

2 16 

3 

1 

1 

4 

281474976710656 

2 48 

5 

18448995933652254721 

4295229439 2 

6 

18013780039499776 

2 16 7 2 74897 2 

7 

7922332684788105606123945 

9841 

281466386710529 2 

8 

5194832314440011219064571 

543158784 

2 60 7 4 23 2 59561 2 

9 

5775028020578783636854257 

0774529 

37 2 701 2 292993041329 2 

10 

2229666183092939997026658 

7262959037621272576 

2 16 19 4 3449 4 4295229439 2 

11 

1463330673647120201450844 

900178197550156472647681 

32363 2 7282397 2 5132726390881 2 

12 

9132817257932632786870155 

6304335790376407269376 

2 48 7 2 i 3 2 ioi77 2 74897 2 259691 2 

13 

6274228310768040852924579 

1977173633010223354344560 

089620481 

3 18 3769 2 15053 2 27205307 2 2607270173 2 

14 

4124060734576866740544330 

9242707472643430878237415 

8659336863744 

2 16 13 2 794009 2 27304061 2 281466386710529 2 

15 

9806175554643239147044270 

0484791252942607177394577 

588251525121 

17489 2 4295229439 2 131825214490835791 2 

16 

1765121615339370515604475 

3103664126594001046686372 

42611924881667927310336 

2 72 7 s 23 2 59561 2 20394769 2 288208447 2 
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